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Abstract 

The constraints on the two Higgs doublet model from the new experimental 
bounds of 5 —> S7 by CLEO and the latest published value of the top quark 
mass by CDF and DO are reanalyzed with the effective Lagrangian covering 
the full QCD corrections from the energy scale of top quark to that of bot¬ 
tom. The reanalysis result shows that the constraints become more stringent 
than that of the earlier analysis, i.e. a bigger region of the parameter space 
of the model is ruled out. 
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It is known that the experimental bounds of 6 ^ S 7 set very strong constraints 
on the two Higgs doublet model (2HDM), a minimal extension of the Standard 
Model (SM). In additional to searching for the neutral Higgs of minimal SM, phe¬ 
nomenologically to investigate possible extensions of SM is also another hot topic in 
particle physics, thus to apply the latest experimental results of the measurement 
on 6 — S 7 |l| and the newly discovery of top quark P] to reexamine the constraints 
on the 2HDM so as to upgrade the allowed values of the model parameters is no 
doubt always to be interesting. 

Reviewing the earlier analysis^], one would hnd that the QCD correction effects 
owing to the change of the energy scale from top quark’s down to that of W boson 
were ignored. Indeed this piece of QCD correction itself is not great, but we treat 
them seriously, and hnally hnd it being not negnigible since this correction affects 
the constraints on the two Higgs doublet model sizable in the report. 

Let us hrst recall the necessary formulas here for later use (mainly from ref. p). 
There are two models for 2HDM to avoid tree-level havor changing neutral currents 
(FCNCs). The hrst (Model I) is to allow only one of the two Higgs doublets to 
couple to both types, u-type and d-type, of quarks [^j but the other doublet is totally 
forbidden by certain discrete symmetry. The second (Model H) is to arrange as that 
one Higgs doublet couples to u-type quarks while the other couples to d-type quarks 
respectively due to a diherent discrete symmetry[|]. It is of interest to note that the 
Model H, as a natural feature, occurs in such a theory as that with supersymmetry 
or with a Peccei-Quinn type of symmetry. 

The piece of the 2HDM Lagrangian for the charged Higgs to the quarks is 
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where V represents the 3x3 unitary Cabibbo-Kobayashi-Maskawa (CKM) matrix. 
Mu and Mu denote the diagonalized quark mass matrices, the subscript L and R 
denote left-handed and right-handed quarks, respectively. For Model I, ^ = V 2 /V 1 ; 
while for Model H, ^ = —vi/v 2 . And ui, V 2 are the magnitude of the vacuum 
expectation values of two Higgs doublets, respectively. 

The relevant effective Hamiltonian after integrating out the heavy top freedom 
is: 

Heff = 2V2GFVthVZY.Ci{l^)Oiif^)- ( 2 ) 

i 

The coefficients Ci{mt) can be calculated from matching conditions at p = m^, and 
Ci{^) can be obtained from their renormalization group equation(RGE): 

= E(7")«C,(,*), (3) 
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where ( 7 ’^)*^- is the anomalous dimension matrix of the operators 0*. If integrating 
out the W boson freedom further, once more six relevant four-quark operators will 
be added 0. We do not repeat all the operators here but only write down the most 
relevant operators: 


Or = {e/ 16 Ti‘^)mbSLa^''bRF^^, 

Os = {g/167r^)mbSLa>^''T%RG%. 


The effective Hamiltonian appears just below the W-scale as 
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For completeness, the explicit expressions of the coefficient of operators at /i = 
Mw are given |Q, 


Gi{Mw) = 0, /or i = 1,3,4, 5, 6; G 2 {Mw) = 1, 
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and together with the coefficients of operators at fi = rrit, 
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= vl/vl, Model J, 
e' = -1, Model II. 

When ^ = 0, the above result (||) reduces to that of SM case0]. 

The running of the coefficients of operators from p = Mw to fi = mb was well 
described in refs.p, |TT|, |T2|, |^. With the running due to QCD, the coefficients 
of the operators at /i = mf, scale: 


C^/f{mb) = + C^iM^) E (9) 

'J i=i 


Here r] = as{Mw)/as{mb), 
hi = 
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Following refs. H, ||, [^, |^, |T3 


BR{B Xs-f)/BR{B X^eV) ~ r(6 ^ s'^)/T{h ceV). 
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Then 


BR{B X,j) ^ \Vt*sVtb\‘^ QaQED N|2 
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where the phase space factor g{z) is given by: 
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g{z) = l-8z^ + 8 / - / - 24/ log 
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here we use rric/mb = 0.316. Afterwards one obtains the B —>■ decay rate 

normalized to the quite well established semileptonic decay rate Br{B Xcev). If 
we take experimental result BR{B XcCV) = 10.8% |T^, the branching ratios of 
B —> Xgj is found to be: 


BR{B Xs'j) ~ 10.8% X 


\V*MHaQED \Cf^{mb)\^ 
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Note that in the above equations the top mass m* is kept as a parameter precisely 
thus we may apply them to computing all the values with various experimental rrit 
as one needs. 

To emphasize the consequences have been ignored in literature, first of all we plot 
the coefficients of the most relevant operators O 7 and Os si fi = Mw versus tan/? 
with and without the QCD corrections of the energy scale running from rrit (with 
the latest experimental value) to Mw in Fig.l for Model I and Model II respectively. 
Here mi, = 4.8GeV, Mw = 80.22GeV and the QGD coupling as{mz) = 0.125 are 
takenUldj. Owing to the mixing of all the relevant operators being small, one may 
see the fact that the effects of the QGD corrections are roughly within ten percent 
and not depend on tan/3 very much. However, one will see soon that the effects, 
though only in ten percent, will make substantial changes for the constraints on 
the parameter space of 2HDM. In order to see the experimental uncertainties of the 
measurements of the branching ratio BR{b —> sy), the top mass and the strong 
coupling constant as{M‘^) how to affect the conclusions, we plot the branching ratio 
BR{h —>■ sy) (the central value also with three standard deviation upper and lower 
bounds achieved by GLEO recently) versus tan/? so as to show the constraints on 
tan/? and on the mass of charged Higgs for Model I (Fig.2a) and Model H 
(Fig.2b), and plot the bends in the parameter space tan/? versus to show the 
constraints caused by taking m* = VJQGeV with ag^M"^) = 0.12 and 0.13 (the 
dashed lines) and as{Mz) = 0.125 with = 163GeH and 189GeH (the solid lines) 
for Model I (Fig.3a) and Model H (Fig.3b) respectively. Here it is interesting to 
note that in Fig.3b, the dashed lines approach to 350GeH and 390GeH respectively 
whereas the solid lines approach to 310GeH and 430GeH respectively when tan/3 
approaches to inhnity. All lines in Fig.3, are obtained by taking the 95% G.L. value 
of GLEO, 1.0x10“^ < Br{B —>• X^y) < 4.2x 10“'^[||]. It is shown that the parameter 
space is more sensitive for changing of mt than for ag, especially in Model H. 

One may see from Figs.2a,3a that for Model I, there are two bands in the tan/?- 
plane, excluded by our reanalysis with the latest measurements on 6 —> s + y 
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and rrit. Finally, as another resnlt of the reanalysis, in Fig.4 we plot the tan/? 
versns plane to show the exclnded region of the parameters for Model II at 
95% C.L.. The solid line corresponds to the conclnsion of the present analysis 
with the central valnes of all the parameters (m^, and the branching ratio of 
the semileptonic decay of b qnark etc.); the dashed line corresponds to a possible 
parameter: V/fep = 0.99 instead of V/fep = 0.950; and the dot- 

dashed line corresponds to the resnlts obtained by other anthorsP|, which do not 
inclnde the piece of the QCD corrections from to Mw, 

In conclnsion, dne to the QCD corrections from rrit to Myy, the new experimental 
valne of rrit and the bonnds for b sy, the constraints for 2HDM are strained 
snbstantially. For instance, the lower bonnd for the mass of the charged Higgs is 
pnt np at least 150GeH for Model II. 
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Figure Captions 


1. Coefficients of operators O 7 , Os, as functions of tan/? in Model I (a) and 
Model II (b). Solid and dashed lines correspond to C 7 with and without QCD 
running from nit to Mw', dash-dotted and dotted lines correspond to Cs with 
and without QCD running from nit to Mw- 

2. Branching ratios ofb^s'y depicted as functions of tan (3, with different masses 
of the charged Higgs in Model I (a), and Model II (b). The CLEO experi¬ 
ment’s central value, upper and lower limit of 95% C.L. for this decay are also 
depicted as solid lines. 

3. tau/d-M^ plane to show excluded region of parameters of Model I (a), Model 
II (b). Every two solid lines denote the allowed and excluded region(s) with 
oisijnz) = 0.125 and nitop = 189GeI/ or mtop = 163GeH; Every two dashed 
lines correspond to that of mMTQGeV and aginiz) = 0.13 or aginiz) = 0.12. 

4. tau/d-M^ plane to show excluded region of parameters of Model II. Solid 
lines corresponds to all central value of parameters; dashed line correspond 
to |I^s^feP/|I46p = 0.99; and dot-dashed line correspond to results without 
QCD Corrections from m* to Mw, obtained by other authors. 
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Fig. 3 (b) 







